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Abstract 

The dynamics of a thin spherically symmetric shell of zero-rest-mass mat- 
ter in its own gravitational field is studied. A form of action principle is used 
that enables the reformulation of the dynamics as motion on a fixed back- 
ground manifold. A self-adjoint extension of the Hamiltonian is obtained via 
the group quantization method. Operators of position and of direction of mo- 
tion are constructed. The shell is shown to avoid the singularity, to bounce 
and to re-expand to that asymptotic region from which it contracted; the 
dynamics is, therefore, truly unitary. If a wave packet is sufficiently narrow 
and/or energetic then an essential part of it can be concentrated under its 
Schwarzschild radius near the bounce point but no black hole forms. The 
quantum Schwarzschild horizon is a linear combination of a black and white 
hole apparent horizons rather than an event horizon. 
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1 Introduction 



According to general relativity, all parts of a massive object definitely disappear if 
the object falls through its Schwarzschild radius. The problem to be tackled in the 
present paper is whether also a quantum system is, or is not, irretrievably lost if it 
falls under its Schwarzschild radius. 

We limit ourselves to a sufficiently simple model so that no approximations are 
needed and the quantum theory can be constructed without problems. In this way, 
an important question about the validity of approximative methods such as WKB 
expansion can also be touched. The simplest system that can ever be invented for 
these aims seems to be a thin shell with its own gravitational field made of light- 
like material, everything spherically symmetric. A Hamiltonian action principle 
|l[ for this system has been transformed to a form suitable for quantization in 
(foregoing paper); this will be used as a starting point. Most of the results of the 
present paper have already been published in a short review 0; here, all derivations 
and calculations will be described in sufficient detail, some new results will be added, 
and a new interpretation of the results will be given. 

The plan of the paper is as follows. In Sec. 2, the starting assumptions and equa- 
tions are collected. The action for the system from Ref. |]J is written down because 
we shall need the form of the constraints. The same action after the transforma- 
tion to a set of embedding variables and Dirac observables is then given, the key 
notion of background manifold is introduced, and the meaning of the new variables 
is discussed. A construction of quantum mechanics including the position and the 
direction-of-motion operators is contained in Sec. 3. The so-called group-theoretical 
quantization method is used, which is well adapted to the problems such as the 
limited ranges of spectra and the construction of a unitary dynamics. The quantum 
mechanics is formulated as a dynamics of the shell on the background manifold; 
this enables straightforward and unique interpretations. In Sec. 4, motion of wave 
packets is investigated. It turns out that no shell reaches the zero radius if it starts 
away from it and so the singularity is avoided. The wave packets contract, bounce 
and then expand, reaching the asymptotic region from which they have been sent in, 
so the dynamics is unitary from the point of view of one family of observers. Some 
of the packets can be sufficiently concentrated near their bouncing point so that an 
essential part of them comes under the corresponding Schwarzschild radius, but no 
event horizon forms. 

In Sec. 5, we consider the seemingly contradictory claims that the quantum shell 
can cross its Schwarzschild radius and still re-expand. The solution of the paradox 
is that if the matter creates a Schwarzschild (apparent) horizon outside then the 
horizon can be, even in the classical version of the theory, of two types: white or 
black, that is, corresponding to the white or black hole horizon in the Schwarzschild 
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spacetime. The "colour" of the apparent horizon in a Cauchy surface depends on the 
direction of motion of the shell: the horizon is black if the shell is contracting and 
it is white if the shell is expanding. The quantum horizon is a linear combination of 
both because the motion of the shell is. The quantum horizon is "grey" , changing 
from mostly black to mostly white. 

The semi-classical approximation fails blatantly near the bouncing point of the 
quantum shell because every classical shell reaches its Schwarzschild radius, forms 
a black hole and falls into the singularity. A cautious discussion of this point is 
given in Sec. 6. In particular, the reason is explained why our results do not prevent 
massive quantum systems from collapsing to black-hole-like objects. 



2 Canonical formalism 

In this section, we shall summarize the formulae derived in Refs. |l| (abbreviated as 
LWF further on) and pf that are needed to start the present paper. 

In LWF, the spherically symmetric metric outside the shell is written in the form 

ds 2 = -N 2 dr 2 + A 2 {dp + N r dr) 2 + R 2 dQ 2 , 

and the shell is described by its radial coordinate p = r. The LWF action reads 

//•OO 
dr pr + y dp(P A A + P R R-H ) , 

and the LWF Hamiltonian is 

H = NH + N^Hp + iVooMoo, 
where := limp^oo N(p), is the ADM energy, H and H p are the constraints, 

AP 2 P\Pr RR" RR'A' R' 2 A 77p 

U = -2R- — + -A A^ + 2A-2 + T 5(/? - r) ' (1) 

Hp = P R R'-P^A-p5(p-r); (2) 

the prime denotes the derivative with respect to p and the dot that with respect to 
r. 

In Ref. 0, the variables rj, r, p, A, Pa, R and Pr have been transformed to the 
embedding variables U(p) and V(p), their canonical conjugates Pu(p) an d Py{p), 
and the shell variables u, v, p u and p v . The pair (u{p), V{p) j defines an embedding 
of the half-axis into the so-called background manifold M. that is covered by the 
coordinates U and V with the ranges 

U + V . n U ~\~V , , 
— - — G (—00,00), G (0,oo). 
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The transformation to embedding variables is determined by a gauge condition, and 
there has been a definite condition used in Ref. 0, where all details are given. 
The background manifold carries then a set of metrics, one representative for each 
geometry. The variables u and v are the coordinates of the shell trajectory in the 
background manifold: 

U = u{t), V = v(t). 
The full action that results has the form of the so-called Kuchaf decomposition 

/f /*oo 
dr (p u u + p v v - npup v ) + Jdrj dp(PuU + P V V - H), (3) 

where H = N u Pu + N v Py; N u (p) and N v (p) are Lagrange multipliers. 

The variables u, v, p u and p v span an extended phase space of the shell. They 
contain all true degrees of freedom of the system. The phase space has non-trivial 
boundaries: 

— it + v 

p u <0, p v <0, >0. (4) 

The constraint surface of the extended action of the shell consists of two components: 
outgoing shells for p v = and in-going shells for p u = 0. 



3 Group quantization 

To quantize the system defined by the action (|3|), we apply the so-called group- 
theoretical quantization method ||. There are three reasons for this choice. First, 
the method as modified for the generally covariant systems by Rovelli || (see also 
and [0]) is based on the algebra of Dirac observables of the system; dependent 
degrees of freedom don't influence the definition of Hilbert space. Second, the group 
method has, in fact, been invented to cope with restrictions such as Eq. (|j). Finally, 
the method automatically leads to self-adjoint operators representing all observables. 

In particular, a unique self-adjoint extension of the Hamiltonian is obtained in 
this way, and this is the reason that the dynamics is unitary. The uniqueness of the 
self-adjoint extension of the Hamiltonian is truly a result of the group quantization 
in the sense that the Hamiltonian operator itself, as calculated from the constraint, 
possesses a one-dimensional family of such extensions. 

To begin with, we have to find a complete system of Dirac observables. Let us 
choose the functions p u , p v , D u := up u and D v := vp v . Observe that u alone is 
constant only along outgoing shell trajectories (p u ^ 0), and v only along in-going 
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ones (p v 7^ 0), but up u and vp v are always constant. The only non vanishing Poisson 
brackets are 

{D u ,p u } = p u , {D v ,p v }=p v . 

This Lie algebra generates a group Go of symplectic transformations of the phase 
space that preserve the boundaries p u = and p v = 0. G is the Cartesian product 
of two copies of the two-dimensional affine group A. 

The group A generated by p u and D u has three irreducible unitary representa- 
tions. In the first one, the spectrum of the operator p u is [0, oo), in the second, p u is 
the zero operator, and in the third, the spectrum is (— oo, 0], see Ref. ||. Thus, we 
must choose the third representation; this can be described as follows (details are 
given in Ref. ||). 

The Hilbert space is constructed from complex functions ip u (p) of p G [0, oo); the 
scalar product is defined by 

f°° dp 

(ip u ,(f) u ) := / — ipM&uip), 
Jo V 

and the action of the generators p u and D u on smooth functions is 



(Pu4>u)(p) = -P4>M, (D u ip u )(p) = -ip- 



dp 



Similarly, the group generated by p v and D v is represented on functions if) v (p)\ the 
group G can, therefore, be represented on pairs [ip u (p),' l Pv(p)) of functions: 



Pu\ipu{p)^v{p)j = (-pif>M,o\ 
Pv(i>u(p),ipv(p)) = ( , -p&v (p) 



d u (Mp),Mp)) = (-v^j^o 
d v {Mp),A(p)) = (°.-^^^)- 

This choice guarantees that the Casimir operator p u p v is the zero operator on this 
Hilbert space, and so the constraint is satisfied. 

Handling the last inequality (ffl) is facilitated by the canonical transformation: 

t = (u + v)/2, r = (-u + v)/2, (5) 

Pt=Pu+Pv, Pr = -Pu+Pv (6) 

The constraint function then becomes p u p v = [p\ — 
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The positivity of r is simply due to its role as the radius of the shell: it is defined 
as a square root of a sum of squares of coordinates with the range R 3 . This suggests 
the following trick. Let us extend the phase space so that r e (—00, +00) and let us 
define a symplectic map / on this extended space by I(t,r,p t ,p r ) = (t, —r,p t , —p r )- 
The quotient of the extended space by / is isomorphic to the original space, and we 
adopt it as our phase space. 

Clearly, only those functions on the extended space that are invariant with respect 
to / will define functions on the quotient. Dirac observables of this kind are, eg., 
p t , pi, the "dilation" D := tp t + rp r = up u + vp v and the square of the "boost" 
J 2 := (tp r + rp t ) 2 = (—up u + vp v ) 2 . The "action" of the map / on the functions p U) 
p v , D u and D v is: 

IpJ = p v , IDJ = D v , IpJ = p u , IDJ = D u . 

There are only two choices for I that preserve these relations in the quantum theory: 

l{^M,^v(p)j = (±^«(p),±^u(p))- 

We choose the plus sign; it is easy to see that the other choice leads to an equivalent 
theory. Observe that the resulting representation of the group G := G <E> (id, /) is 
irreducible. 

There are two eigenspaces of /: one to the eigenvalue +1, consisting of the pairs 
with i]) u (p) = ip v {p)i the other to the eigenvalue —1, containing the pairs with 
i>u(p) — ^"^vip)- If we choose one of these eigenspaces as our final Hilbert space, 
we obtain a representation of the classical algebra on the quotient space. Again, the 
two possible choices give equivalent theories. The final result can easily be brought 
to the following form. The states are determined by complex functions <p(p) on R + ; 
the scalar product (</?, ip) is 

fOO 



Jo V 



let us denote the corresponding Hilbert space by /C. The representatives of the above 
algebra are 

(Pt<p)(p) = -P<f(p), 

(Pr^)(p) = PV(P)» 

The next question is that of time evolution. Time evolution of a generally co- 
variant system described by Dirac observables may seem self-contradictory or gauge 
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dependent. Here, we apply the approach that has been worked out in P and || 
using the symmetry group of time shifts found in Sec. 2 of Ref . |2| , which is gener- 
ated by the function p t . The operator — p t has the meaning of the total energy M of 
the system. We observe that it is a self-adjoint operator with a positive spectrum 
and that it is diagonal in our representation. The parameter t of the unitary group 
U{t) that is generated by — p t is easy to interpret: t represents the quantity that is 
conjugated to p t in the classical theory and this is given by Eq. (|5|). Hence, U(t) de- 
scribes the evolution of the shell states between the levels of the function {U + V)/2 
on M. 

The missing piece of information of where the shell is on Ai is carried by the 
quantity r of Eq. (Q). We try to define the corresponding position operator in three 
steps. 

First, we observe that r itself is not a Dirac observable, but the boost J is, and 
that the value of J at the surface t = coincides with rpt- It follows that the 
meaning of the Dirac observable Jp^ 1 is the position at the time t = 0. This is in a 
nice correspondence with the Newton- Wigner construction on one hand, and with 
the so-called evolving constants of motion by Rovelli |7J on the other. 



Second, we try to make Jp^ 1 into a symmetric operator on our Hilbert space. As 
it is odd with respect to /, we have to square it. Let us then chose the following 
factor ordering: 

r 2 := — J- J- = -VP—-. (7) 

Vp p Vp d P 2 Vp 

Other choices are possible; the above one makes f 2 essentially a Laplacian and this 
simplifies the subsequent mathematics. Indeed, we can map /C unitarily to L 2 (R + ) 
by sending each function ip{p) G /C to ip{j>) G L 2 (R + ) as follows: 

V>(p) = -pV(p)- 
VP 

Then, the operator of squared position f 2 on L 2 (R + ) corresponding to f 2 is 

f 2 = -^r 2 (^( P )) = = ~M(p). 

Third, we have to extend the operator f 2 to a self-adjoint one. The Laplacian 
on the half-axis possesses a one- dimensional family of such extensions ||rTJ. The 
parameter is a G [0, 7r) and the domain of A Q is defined by the boundary condition 
at zero: 

?/>(0) sin a + ^'(0) cos a = 0. 
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The complete system of normalized eigenf unctions of A a is given by: 
~ / 2 r cos a cos rp — sin a sin rp 

M r ,p) = \ -, = „ . =2= — ; 

if a G (0, 7r/2), there is one additional bound state, 

$a(b,p) = 7= exp(-ptana), 
V 2 tan a 

so that 

-A a 4i a (r,p) = r 2 ip a (r,p), 
-A a ^ a (b,p) = -tan 2 a ip a (r,p). 

The corresponding eigenfunctions ip a of the operator f 2 are: 

' 2p r cos a cos — sin a sin rp 



^ a (r,p) 




71 V r 2 cos 2 a + sin 2 a 



and we restrict ourselves to a G [vr/2, 7r], so that there are no bound states and the 
operator f is self-adjoint. 

To restrict the choice, we apply the idea of Newton and Wigner. First, the 
subgroup of G that preserves the surface t = is to be found. This is, in our 
case, C/d(A) generated by the dilatation D. Then, in the quantum theory, the 
eigenfunctions of the position at t — are to transform properly under this group; 
this means that the eigenfunction for the eigenvalue r is to be transformed to that 
for the eigenvalue Ur,(X)r, for each A. The dilatation group generated by D acts on 
a wave function ip{p) as follows: 

^(p)^U D (\)4>(p)=4>(e~ x p), 

where £/d(A) is an element of the group parameterized by A. Applying £/d(A) to 
*Pa(r,p) yields 

, . . x/9 j2p r cos a cos(e~ x rp) — sin a sin(e _A rp) 

U D (X)Mr,p) = e- A/2 \/-r 7= = , . 2 — -■ 

V 7T cos 2 a + sm a 

The factor e _A / 2 in the resulting functions of p keeps the system ^-normalized. 
Let a = 7r/2; then 

U D (\)^ /2 (r,p) = e- x/2 ^ /2 (e~\p). 

Similarly, for a = it, 

U D (X)^(r,p) = e- A/2 ^(e" A r,p), 
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but such relation can hold for no other a from the interval [tt/2, 7r], because of the 
form of the eigenfunction dependence on r. Now, Newton and Wigner require that 



U D (X)^(r,p) = e- x ^( 



e r,p). 



Then all values of a except for a = tt/2 and a = ir are excluded. 

We have, therefore, only two choices for the self-adjoint extension of f 2 : 



The construction contains a lot of choice: the large factor-ordering freedom, and 
the freedom of choosing the self-adjoint extension. One can react to this ambiguity 
in two ways. 

The first is to ask how the different choices influence the results. It seems plausible 
that the qualitative, rough properties of the quantum system will be the same for 
all possible choices. We hope (provisionally) that this is true. 

The second question to ask is how the position is, in fact, measured in praxis. 
This question hits the crux of the problem. Indeed, the Newton- Wigner construction 
may be formally elegant but, to my knowledge, nobody managed to describe the 
corresponding measurement. If we search for methods of how the position of various 
constituents in a microscopic system is measured, we find the scattering method to 
dominate. For that it is necessary to use a particular coupling the system under 
study, a crystal, say, has with another agent, X-rays, say. One has to send the 
X-rays onto the crystal and to view what comes out. 

It seems, therefore, that the following approach would be more reliable than 
attempts at a formal definition of a position operator of the shell. One can try, for 
example, to couple the shell to some field, the quanta of which could be emitted by 
the shell on its way down and up. The quanta will, or will not reach the asymptotic 
observers and their properties at infinity might reveal something of what is going 
on with the shell. This is a future project because it will be mathematically more 
difficult than our provisional attempt with the position operator. 

Another observable that we shall need is fj; this is to tell us the direction of motion 
of the shell at the time zero, having the eigenvalues +1 for all purely outgoing shell 
states, and —1 for the in-going ones. In fact, in the classical theory, 77 = — sgnp r , 




r > 0, 




and 




Let us select the first set, Eq. (||) 
is finished. 



by that, the construction of a position operator 
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but p r does not act as an operator on the Hilbert space /C, only p 2 T . Hence, we need 
the following trick. 

Consider the classical dilatation generator D = tp t + rp r . It is a Dirac observable; 
at t — 0, its value is rp r . Thus, for positive r, the sign of —D at t — has the 
required value. On the quotient space, the values at negative r correspond to the 
/-mapped states with positive r, and, as D is /-invariant, the relation of the sign to 
the direction of motion is again valid. Hence, we have the relation: 

sgnZ) = -77 t=0 . 

The normalized eigenfunctions ip a (p) of the operator D are solutions of the dif- 
ferential equation: 

D4>a{p) = aip a (p). 

The corresponding normalized system is given by 



Hence, the kernels P±(p,p') of the projectors P± on the purely out- or in-going states 
are: 

P + M)= f daMp)^f-, P-(P,P)= r daUP)^ 1 

J -oo P JO P 

so that 

P oo 

(#)(p) = / dp'[p + (p,p') - p_(p,p')]^(p')- 



This finishes our construction of the shell quantum mechanics. 



4 Motion of wave packets 

We shall work with the family of wave packets on the energy half-axis that are 
defined by 

where k is a positive integer and A is a positive number with dimension of length. 
Using the formula 



'o 



dpp n e-^ = —, (9) 
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which is valid for all non-negative integers n and for all complex v that have a 
positive real part, we easily show that the wave packets are normalized, 



1. 



f°°dp 
Jo P 
The expected energy, 

m kx -.= r^», 

Jo P 

of the packet can be calculated by the same formula with the simple result 

The (energy) width of the packet can be represented by the mean quadratic devia- 
tion, AM K \, which is 



Hence, by choosing k and A suitably, we can approximate any required energy and 
width arbitrarily closely. 

The time evolution of the packet is generated by —pt- 

Let us calculate the corresponding wave function ty K \(r, t) in the r-representation, 

dp 



Vnxfar) := / 
Jo 



P 



where the functions tp(r,p) are defined by Eq. (pj). Formula @ then yields 

i K\(2\y +l ' 2 



It follows immediately that 



(\ + it + ir) K+1 (\ + it-ir) K+1 



(10) 



lim|# KA (t,0)| 2 = 0. 

! H) 

The scalar product measure for the r-representation is just dr because the eigen- 
functions @ are normalized, so the probability to find the shell between r and r + dr 
is \^ KX {t } r)\ 2 dr. 

Our first important result is, therefore, that the wave packets start away from the 
center r = and then are keeping away from it during the whole evolution. This 
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can be interpreted as the absence of singularity in the quantum theory: no part of 
the packet is squeezed up to a point, unlike the shell in the classical theory. 

Observe that the equation ty K \(t,0) = is not a result of a boundary condition 
imposed on the wave function. It is a result of the unitary dynamics. The nature 
of the question that we are studying requires that the wave packets start in the 
asymptotic region so that their wave function vanishes at r = for t — ► — oo; this 
is the only condition put in by hand. The fact that the dynamics preserves this 
equation is the property of the unique self-adjoint extension of the Hamiltonian 
operator. 

A more tedious calculation is needed to obtain the time dependence f K \(t) of the 
expected radius of the shell, 



r K x(t) := / drr\V KX (t,r)\ 2 . (11) 
J o 

Let first k = . The wave function of the packet then is 

* llA (f.r) = 2, /A 



vr r 2 + (A + it) 2 ' 

so the expectation value of f is 

M 4Af, r 3 
ro\(t) = — / dr 



vr J (r 2 + \ 2 -t 2 ) 2 + 4\H 2 ' 
This integral diverges logarithmically, so 

f x(t) = oo. 

Let The substitution of Eq. (|l(J) into ( |TT|) leads to: 

where 

x f 1 1 



4a (*) = / rdr 
Jo 

/•OO 

4a (*) = / rdr 



[(r + t) 2 + A 2 ]* +1 [(r - t) 2 + A 2 ] K + X 
1 1 



[(A-ir) 2 + t 2 ] K + 1 [(A + ir) 2 + t 2 ] K + 1 
The first integral can be brought by elementary methods to the following form: 

T / \ i i y* 

I ^(t) = - 7 - r —^- + t 



K{t 2 + \ 2 Y J- t (s 2 + A 2 ) K + r 
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Let us calculate the second integral. We obtain after a simple rearrangement: 

r + iX iA 



•1 



+ 



(If 

r — iX 



o 



[(r + iA) 2 -t 2 ]^ 1 [(r + zA) 2 -t 2 ] K+1 
iA 



[(r - iA) 2 - t 2 ] K+1 [(r - iA) 2 - t 2 ] K+1 J ' 



This suggests the introduction of integration contours C\ defined in the complex 
plane by z = r + iA for r G (0, oo), and C2 by z = r — iA, r e (0, 00). Then J K \(t) 
can be written as follows: 



Ci 



c 2 



iA 



^2_ f 2) K +l (^2_ t 2)*+l 

z iA 
+ 



The integrals of the first terms in the square brackets can be done immediately: 



J K x(t) 



1 



1 



\K+l 



K(x 2 + t 2 y 

We obtain as the final result: 



iA 



dz 



-C!+C 2 



> 2 _ t 2) K+ l 



:i2) 



f KX {t) 



1 (d) 2 (2A) 



2k+1 



1 



+ t 



dx 



2tt (2k)! [K{X 2 + t 2 Y 

iX{ I)"' 1 



+ zA(-l) K+1 / 
Jc 



t (x 2 + A 2 ) K + X 
dz 



Cl yz 2 - t 2 )^ 



C 2 



( z 2 _ py+i 



(13) 



In fact, the R. H. side diverges for k = so, in this sense, this formula can be 
considered as completely general, ie., valid for all k and t. 

Let us study some properties of the function f K \(t). Eq. (|13|) implies that 

f K \(t) = f K \(-t), 

so the average motion of the packet is symmetric under time reversal. Eq. ( JTB| ) 
is also suitable for the calculation of the expansions about the points t = and 
t = ±00. Consider first the point t = 0. Expanding the first term in the square 
bracket is easy: 



n(X 2 + t 2 Y kX 2 



2k 



The series on the R. H. side converges for \t\ < X. 
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To expand the next term, we expand the integrand in the powers of x/X; the 
series converges for \t\ < X. Integrating term by term yields: 



dx 



t (x 2 + A 2 ) K+1 A : 



k=0 



K + k 

2k + 1 V k 



2k+2 



Again, this series converges for \t\ < X. 

A similar method can be applied to the remaining integrals: 

y2K+2 ( 

fc=0 V 



1 

(z 2 - t 2 ) K+1 ~ ~ 2 



k + k 
k 



The convergence is granted for |t| < \z\. As the minimal \z\ along both contours is 
A, the expansion is always valid for |t| < A. Then 

ix(-ir +i ' 



\K+1 



\2k L^i 



(-1)' 



X 2k ^2n + 2k + l 

k=0 



Collecting all terms, we obtain the expansion around t = 0, 



k + k 
k 



r K x{t) 



X (k\) 2 2 



2o2k+1 



7T (2k)\ 



K + l 



+ 



k(2k+ 1) 
(-l) k+1 



^ {2k + 2k + l)Jfe(2ife - I] 



K + k — 1 
k 



2k 



(14) 



and the equation holds for \t\ < X. As the k = 1 term in Eq. (ITJJ) is positive, there 
is a minimal expected radius f K \(0) at t = 0, 



fnx(0) 



12 IK { k\Y k + 1 X 



> 0. 



(15) 



7T {2k)\ k k + 1/2 

For a large k, the minimum at t = may be only local and/or the curve may oscillate 
for t G (-A, A). 

Let us turn to the asymptotics t — > ±00. It is sufficient to consider the case 
i — » 00 because the other one is obtained by t 1— > — t. The first term in the square 
brackets in Eq. ( |13| ) is clearly of the order 0(t _2 ' t ) and it is not difficult to convince 
oneself that the last two terms are both of the order 0(t~ 2K ~ l ). 

The second term need more care. First, we use the relation 



dx 



(x 2 + X 2 )^ 1 



so that we obtain 



dx 



(x 2 + X 2 ) K+1 k\ 



1 d 
2XdX 



1 d 



K r t 



dx 



x 2 + A 2 



2A dX 



— - I — — : — 7— — arctan— 



t 



t 



X 



X 



13 



For t/X > 0, the following formula holds: 

t IT X 

arctan— = arctan— . 

A 2 t 

Using it and expanding the function arctan(A/t) around zero leads to 

t r dx = — f-i^yi- 2 r-iv ( d vf^ 2k 

J_ t (x 2 + X 2 )-+ 1 2 k k\\ XdXj X 1 ' \d{Xy) f^2k + l\t 

The series converges for t > X and can be differentiated term by term in this interval. 
The first non zero term comes only from k = k and it has the value 

2k! 1 



It holds also 



2k + 1 t 2K 
1 d Y 1 (2k- 1)!! 



XdX X X 2K+1 



Hence, 



t dx _ (2k -1)!! _ 2 
J_ t (x 2 + X 2 ) K+1 ~ irt 2"\*"+ 1 K\ + [t } 



Substituting this into Eq. (|13|) and using the symmetry t ^ —t, we obtain for both 
cases t — > ±oo: 

f KX (t)^\t\ + 0(r 2K ). (16) 

A further interesting question about the motion of the packets is about the portion 
of a given packet that moves in — is purely in-going — at a given time t. The portion 
is given by ||-P_-?/v\|| 2 , where P_ is the projector defined in Sec. |. Let us calculate 
this quantity. 

If we write out the projector kernel and make some simple rearrangements in the 
expression of the norm, we obtain: 

/OO /"OO / poo \ 

dq' / dq" / dara(e q, )Me q ") e 9 ")^ A (t, e"'), 

oo J-oo \Jo J 

where the transformation of integration variables p' and p" to e q ' and e q " in the 
projector kernels has been performed. 
The integral in the parenthesis, 

dar a (e q ')Me g ") = 7T / dae*** ~*\ 
2-n Jr. 
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is a kernel in an integral that is exponentially damped at the infinities. Thus, we 
can calculate it limit, 



oo 



lim / dae ia{q "- q,) - ea = — lim 



2ir e^o J Q 2ir e^o (q" - q') + ie 

2n q" — q' 2 

where V denotes the principal value. 

Doing the integral over the ^-function gives the simple result: 

1 . . 1 

The rest can be written as follows: 

\\P-^x\\ 2 = \+ 

i r r uu. .Vv,. ,^cost(g"-g')+isint( g "-g') 



dq'V / dq"iP KX {e q )^ KX {e 



2tt y_ 00 y_ 00 q" - q' 

The integrand is a sum of a symmetric and an anti-symmetric functions of the 
variables q' and q" . The principal value integral annihilates the anti-symmetric 
part. The integral from the symmetric part is already regular, and we can write the 
final formula: 

hp / lis 1 1 Hj'Tj"! (I'm < ^ sint(g"-g') n7 , 
II p -^a|| =9-^-/ d 1 / ^ ^A(e 9 )^ KA (e 9 ) — - — . (17) 

Let us calculate the in-going portion for some simple values of t. Thus, for t — 0, 
we obtain immediately: 

1 

2" 

At the time zero, the probabilities to catch the shell going in or out are equal. 
The limit t — > ±00 can be obtained, if we use the formula: 



Hence, 



sintx 
lim = ±tio[x) 

t->±oc x 



lim Sint(e9 "- eg,) = ±n- -6(e*' - e*). 

t^ioo q" - q> q" - q' 



Substituting this into the integral of Eq. fll7|) and returning back to the variables p' 
and p" results in: 

116/ 112 1 1 f°° d P' f°° d P" , ( >M P"-P' x , n n 

2 2 J p' J p" logp"-logp' 
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The expression 



\ogp" — \ogp' 



is smooth and equal to p' at p" = p' . Hence, finally 



p^x\\L-oo = i, \\p-M\Loo = o, 



and we have only in-going, or only outgoing shells at the infinity. 

The obvious interpretation of these formulae is that quantum shell always bounces 
at the center and re-expands. We can, however, ask further questions. For example, 
what is the time delay of the re-expansion as compared, say, with the same trajectory 
in the background manifold M that carries the flat metric 



in our coordinates U and V? To find this time delay, the "true" metric with respect 
to these coordinates had to be calculated. The metric is determined by the quantum 
state in a similar way as the position and the colour of the horizon are (see the next 
section). However, unlike the points and the metric inside Ai, the points and metric 
in the asymptotic region are gauge invariant quantities. The method by which we 
should calculate the asymptotic metric ought to make the gauge invariance of the 
result transparent. Such a method has first to be developed. 

The result that the quantum shell bounces and re-expands is clearly at variance 
with the classical idea of black hole forming in the collapse and preventing anything 
that falls into it from re-emerging. It is, therefore, natural to ask, if the packet 
is squeezed enough so that an important part of it comes under its Schwarzschild 
radius. We can try to answer this question by comparing the minimal expected 
radius f KA (0) with the expected Schwarzschild radius f K \ H of the wave packet. The 
Schwarzschild radius is given by 



where M P is the Planck energy. Now, the values of k and A for which a large part 
of the packet gets under its Schwarzschild radius clearly satisfy the inequality 



ds 1 = -dUdV + (1/4) (-17 + Vfdtt 2 



(18) 



f KXH = 2GM KX 



M kX 



My 



r K \{0) < r KXH , 



or 



(AM P ) 2 < 2tt 



k(k+ 1/2) 2 (2k)\ 
k+1 2 2K (/d) 2 ' 



(19) 
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Interpreting A roughly as the spatial width of the packet, we have XMp 3> 1 for 
reasonably broad packets. Then the right-hand side can be estimated by the Stirling 
formula: 

kU+1/2) 2 (2k)! , — 
k + 1 2 (re!) 



Substituting this into the inequality ( |19|) yields 



M KX > ^M P , (20) 



/2tt 

which implies that the threshold energy for squeezing the packet under its Schwarz- 
schild radius is much larger than the Planck energy. For narrow wave packets, we 
have that XM P « 1, so the inequality (fL9|) is satisfied, and the threshold energy 
is about one Planck energy. The inequality (p0|) expresses, therefore, always the 
desired property. To summarize: Reasonably narrow packets can, in principle, get 
under their Schwarzschild radius; their energy must be much larger than Planck 
energy. Even in such the shell bounces and re-expands. 

This apparent paradox will be explained in the next section. 



5 Grey horizons 

In this section, we try to explain the apparently contradictory result that the quan- 
tum shell can cross its Schwarzschild radius in both directions. The first possible 
idea that comes to mind is simply to disregard everything that our model says about 
Planck regime. This may be justified, because the model can hardly be considered 
as adequate for this regime. However, the model is mathematically consistent, sim- 
ple and solvable; it must, therefore, provide some mechanism to make the horizon 
leaky. We shall study this mechanism in the hope that it can work in more realistic 
situations, too. 

To begin with, we have to recall that the Schwarzschild radius is the radius of 
a non-diverging null hyper-surface; anything moving to the future can cross such 
a hyper-surface only in one direction. The local geometry is that of an apparent 
horizon. (Whether or not an event horizon forms, that can also depend on the 
geometry near the singularity [0). However, as Einstein's equations are invariant 
under time reversal, there are two types of Schwarzschild radius: that associated with 
a black hole and that associated with a white hole. Let us call these Schwarzschild 
radii themselves black and white. The explanation of the paradox that follows from 
the model is that quantum states can contain a linear combination of black and 
white horizons, and that no event horizon forms. We call such a combination a grey 
horizon. 
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The existence of grey horizons can be shown as follows. The position and the 
"colour" of a Schwarzschild radius outside the shell is determined by the spacetime 
metric. For our model, this metric is a combination of purely gauge and purely 
dependent degrees of freedom, and so it is determined, within the classical version 
of the theory, by the physical degrees of freedom through the constraints. 

To explain the idea in more detail, less us start with the general case in the 
ADM formalism. There are 16 canonical variables, the 6 components of the three- 
metric qki, the 6 components of the conjugate momentum 7i kl , 1 lapse and three 
shift functions. These can be decomposed (non-uniquely) into physical, gauge, and 
dependent variables. Fixing the gauge variables by hand (this also includes some 
boundary conditions in non-compact cases) means that a particular space-like sur- 
face E is chosen, and a particular coordinate system x k , k = 1,2,3, is lain onto 
this surface. Then the constraints turn into differential equations determining the 
dependent part of q k i and n kl in terms of the physical one and so the tensor fields 
qki and 7i H are determined uniquely along E in the coordinates x k by the physical 
degrees of freedom. By this, the full spacetime metric g^ v and all its first derivatives 
are known at each point of E. Indeed, if we choose the Gaussian coordinates x°, x k , 
adapted to E, then the four- metric at E is 

ds 2 = -(dx ) 2 + q kl dx k dx l , 

and the derivatives of this metric with respect to the coordinates x° and x k are given 
by 

dgoo = Q ^oo = Q dgok 

dx° dx k dx° 

dgm = _ 2K dg H = dg kt dg ok 

dx° ~ kh dx m ~ dx m ' dx 1 
where K k i ■— (det q k i)~ 1 ^ 2 (l/2 q mn 7r mn q k i — 7Tjw) is the second fundamental form of 
the surface E. Observe that the choice of Gaussian coordinates is equivalent to 
specifying the lapse and shift at E by hand. The lapse and shift could also be fixed 
by the condition that the gauge is preserved by the evolution. 

Let S be a closed two-surface on E. We can calculate the Gaussian coordinates 
adapted to S in E from the metric q k i of E; let they be x' A and x' 3 , A = 1, 2, so that 
S is given by x' 3 = and x' 3 increases in the outside direction. Let the corresponding 
components of the tensor fields be q kl and K' kl . Then the induced two-metric on S 
is q' AB and the full three-metric on E is 

ds 2 = (dx' 3 ) 2 + q' AB dx' A dx' B . 

Now, let and n M be null vectors orthogonal to S , Z M being the outgoing and n M 
the in-going one. Their component in terms of the coordinates x'° := x° and x' k are 

= (1,0,0,1), n'" = (1,0, 0,-1). 
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Then 



d 9AB m _ 9 9ab , 9 9ab _ 2K > ,^ab 
dx'» dx'° dx' 3 AB dx' 3 ' 

and, similarly, 

dx'» " ab dx , 3 , 

We can, therefore, check, whether or not the following equation holds 

9 '-(-2/^ B± ^f)=0, (21) 

and so can find, if S is an out- (in-)going apparent horizon — Eq. (^l|) is then valid 
with the above (lower) sign. 

Let us look to see how this algorithm works for the shell model. The constraints 
are Pu{p) — and Pv{p) = 0. If the transformation from the original variables 
A (p), R(p), p x(j>), Pr(p), V, r and p to U(p), V(p), P v {p) P v (p), V, u i v i Vu and p v 
were known, it would provide the functionals: 

P u (p)=P u [\,R,P x ,P R , V ,v,p;p), 

and 

P v (p)=P v [X,R,P x ,P R ,r ] ,r,p;p). 

The transformation is not known explicitly, but we know that the constraint equa- 
tions 

Pu[\ R, P A , Pr, V, r, p; p) = 0, P v [A, R, P x , P R , V , r, p; p) = (22) 

are equivalent to the original constraints, Eqs. ([!]) and @). Hence, our first trick is 
to work with Eqs. (JT|) and (H) instead of Eqs. (|22]). 

The constraints (|l|) and (0) contain the physical variables rj, M u and v also 
through r and p. We can choose the gauge variables to be R(p) and A(p). A fixed 
function R(p) determines p in terms of the geometrical quantity R and so it fixes 
a radial coordinate along S. A(p) contains derivatives of the embedding functions, 
so it determines the slope of the embedding at each p. Integrating the slope gives a 
family of surfaces; a suitable boundary condition at infinity selects one of them. 

In order to obtain a suitable surface £ the functions R(p) and A(p) have to 
satisfy some further boundary conditions at the infinity, at the shell and at the 
regular center. The condition at the infinity, p — > oo, is to guarantee that E is 
asymptotically flat. That at the shell is necessary in order that E is smooth across 
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the shell. Finally, at p = 0, we require that E cut the regular center rather that the 
singularity and that it be a smooth surface at this point. Similarly, Pr and Pa have 
to satisfy suitable boundary conditions at p = 0, p = r and p — > oo. The explicit 
form of these boundary conditions are carefully discussed in [|[] . 

Then the constraints become equations for the two functions Pr(p) and P\(p)- 
Eq. (0) is an algebraic equation for P#; solving it and inserting the solution into Eq. 
(0) gives an ordinary differential equation for Pa. The differential equation, together 
with the boundary conditions, determines Pv(p) uniquely, and this, in turn, together 
with Eq. (P, gives Pr(p). The solution is unique. From the known functions R(p), 
A(p), Pr(p) and P\(p), we can determine qui and K k i along E and check Eq. (|2ip . 

One can try to solve Eqs. ([[]) and (0) for Pr(p) and P\(p) explicitly, by choosing 
the functions R(p) and A(p) in some way that simplifies the equations. Instead, we 
use the uniqueness of the solution in the following simple trick. Any solution of the 
constraint equations in the spherically symmetric case defines an initial data and 
surface for a solution to Einstein's equations that is itself spherically symmetric. 
Hence, every such solution of constraints forms a space-like surface that can be 
embedded in some Schwarzschild spacetime. There will always be the Schwarzschild 
solution of mass zero inside the shell, and the Schwarzschild solution of mass M 
outside it. 

In this way, we find by inspection from the Kruskal diagram: If the shell is in- 
going, n = — 1, then it is contracting and any space-like surface containing such a 
shell can at most intersect an outgoing apparent horizon at the radius R = 2GM, 
independently of which of the two infinities the surface is connecting the shell with. 
Analogous result holds for rj = +1, where the shell is expanding. The corresponding 
Ph is determined by the equation R(ph) — 2GM, and the horizon will cut E if 
and only if pu > r. We can assign the value +1 (—1) to the horizon that is out- 
(in-)going and denote the quantity by c (colour: black or white hole). Then c = — rj. 

In particular, if we choose the gauge so that R(p) = p, then r is just r(t), where 
t is the value of the parameter t at which the shell intersect E, and we have: 

1. The condition that an apparent horizon intersects E is r t < 2GM. 

2. The position of the horizon at E is p# = 2GM. 

3. The value of c is c = —rj. 

In this way, questions about the existence and colour of an apparent horizon 
outside the shell are reduced to equations containing dynamical variables of the 
shell. In particular, the result that c = — r] can be expressed by saying that the shell 
always creates a horizon outside that cannot block its motion. All that matters is 
that the shell can bounce at the singularity (which it cannot within the classical 
theory). 
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These results can be carried over to quantum mechanics after quantities such as 
2GM — r and rj are expressed in terms of the operators describing the shell. Then 
we obtain a "quantum horizon" with the "expected radius" 2GM and with the 
"expected colour" —f\ to be mostly black at the time when the expected radius of 
the shell crosses the horizons inwards, neutrally grey at the time of the bounce and 
mostly white when the shell crosses it outwards. 

This proof has, however, two weak points. First, the spacetime metric on the 
background manifold is not a gauge invariant quantity; although all gauge invariant 
geometrical properties can be extracted from it within the classical version of the 
theory, this does not seem to be possible in the quantum theory [[□]. Second, 
calculating the quantum spacetime geometry along hyper- surfaces of a foliation on 
a given background manifold is foliation dependent. For example, one can easily 
imagine two hyper-surfaces E and E' belonging to different foliations, that intersect 
each other at a sphere outside the shell and such that £ intersects the shell in its 
in-going and E' in its outgoing state. Observe that the need for a foliation is only 
due to our insistence on calculating the quantum metric. 

The essence of these problems is the gauge dependence of the results of the 
calculation. However, it seems that this dependence concerns only details such as 
the distribution of different hues of grey along the horizon, not the qualitative fact 
that the horizon exists and changes colour from almost black to almost white. Still, 
a more reliable method to establish the existence and properties of grey horizons 
would require another material system to be coupled to our model; this could probe 
the spacetime geometry around the shell in a gauge-invariant way. 

It may still seem difficult to imagine any spacetime that contains an apparent 
horizon of mixed colours. Nevertheless, examples of such space-times can readily be 
constructed if the assumption of differentiability is abandoned. A continuous, piece- 
wise differentiable spacetime can make sense as a history within the path integral 
method. 

The simplest construction of this kind is based on the existence of the time 
reversal isometry T as defined in the foregoing paper []2| that maps an in-going shell 
spacetime onto an outgoing one. 

Let us choose a space-like hyper-surface Ei crossing the shell before this hits the 
singularity in a (— 1, M, u)-spacetime A4, and find the corresponding surface TEi 
in the spacetime TM with the parameters (1, M, v). Then we cut away the part of 
Ai that lies in the future of Ei and the part of TM in the past of TE^ As their 
boundaries are isometric to each other, the remaining halfs can be stuck together in 
a continuous way. In the resulting spacetime, the shell contracts from the infinity 
until it reaches Ei at the radius r\\ then, it turns its motion abruptly to expand 
towards infinity again. There is no singularity and the spacetime is flat everywhere 
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inside the shell. If r\ < 2GM, then there is an apparent horizon at R = 2GM. 
It comes into being where the in-going shell crosses the radius r = 2GM and is 
outgoing (black) until it reaches E^ Then, it changes its colour abruptly to white 
(in-going) and lasts only until the outgoing shell crosses it again. 

The space-like hyper-surface Si can be chosen arbitrarily in M.. The construc- 
tion can, therefore, be repeated in the future of TEx in an analogous way so that 
we obtain a spacetime with two "pleats"; the shell contracts, then expands, then 
contracts again and hits the singularity. The horizon starts as a black ring, then 
changes to a white one, and then it becomes black for all times. This history is, 
however, not continuous. Clearly, one can repeat the construction arbitrary many 
times; this leads to a "pleated" spacetime with a zig-zag motion of the shell and 
alternating horizon rings of white and black colour. If the spacetime is to be singu- 
larity free, however, there must be an odd number of pleats and an even number of 
rings, beginning with the black ring and ending with a white one. 

The conditions that the surface E x cuts the trajectory of the shell at some small 
value of the Schwarzschild radial coordinate R, is smooth and space-like everywhere 
and hits the space-like infinity i° for large values of R allow a considerable freedom. 
We can require in addition that Ei joins smoothly to the surface T — Ti, where T is 
the Schwarzschild time coordinate and Ti some constant so that Ei coincides with 
T = Ti for all values of R larger than, say, R\. It is clear from the Penrose diagram 
that such a Ex can ran arbitrarily close to the incoming shell trajectory and can be 
joined to T = Ti for arbitrary low value of T\ G (—00, 00), if R\ is chosen sufficiently 
large. On the other hand, for R\ = 2GM + e, Ei can join T — Ti for arbitrarily 
large Ti G (—00, 00), just if e > 0. 

Consider now an observer at the fixed value R of the Schwarzschild radius in 
each shell spacetime. We shall choose Ei in such a way that R\ < R . With this 
choice, the observer trajectory R = R remains smooth at Ei. Then, the lower 
bound on the possible values of Ti is T c , which is the Schwarzschild time of the 
point at which the observer crosses the shell. There is, however, no upper bound 
on Ti. Hence, we can construct a one-pleat spacetime for each value of Ti from the 
interval (T c , 00) with a smooth trajectory of the observer. For each value of Ti, the 
observer will measure the proper time 



between his two encounters with the shell. Thus, the time delay can be made 
arbitrarily small or large. (Of course, all such histories and many others must 
be integrated with some suitable measure in a path integral in order to obtain a 
reasonable value of the delay). 
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Let us choose a gauge in each spacetime constructed above such that the trajec- 
tory of the shell is V = u for the in-going part and U = u for the outgoing one. 
Then, the metric in the asymptotic region, where the observer is, will read 

ds 2 = -A(U, V)dUdV + R 2 (U, V)dR 2 , 

and it is clear that the functions A(U, V) and R(U, V) must have different forms for 
different values of T±, or else the proper time At measured by the observer will be 
independent of T\. In most cases, the asymptotic behaviour of the metric in this 
gauge will be different from (]X8f) . On the other hand, in each such spacetime, there 
will be double null coordinates U\ and V\, say, in which the metric will have the 
asymptotic behaviour fll8). However, the trajectory of the shell will then be given, 



with respect to the coordinates U\ and V\, by different equations for different values 
of T x . 



6 Concluding remarks 

Comparison of the motion of wave packets of Sec. 4 with the classical dynamics 
of the shell as described in Sec. 3 of shows a marked difference. Whereas all 
classical shells cross their Schwarzschild radius and reach the singularity in some 
stage of their evolution, the quantum wave packets never reach the singularity, but 
always bounce and re-expand; few of them manage to cross their Schwarzschild 
radius during their motion. This behaviour is far from being a small perturbation 
around a classical solution if the classical spacetime is considered as a whole. Even 
locally, the semi-classical approximation is not valid near the bouncing point. It is 
surely valid in the whole asymptotic region, where narrow wave packets follow more 
or less the classical trajectories of the shell. 

The most important question, however, concerns the validity of the semi-classical 
approximation near the Schwarzschild radius. We have seen that the geometry near 
the radius can resemble the classical black hole geometry in the neighbourhood of 
the point where the shell is crossing the Schwarzschild radius inwards. Then, the 
radius changes its colour gradually and the geometry becomes very different from 
the classical one. Finally, near the point where the shell crosses the Schwarzschild 
radius outwards, the radius is predominantly white and the quantum geometry can 
be again similar to the classical geometry, this time of a white hole horizon. 

If the change of colour is very slow then the neighbourhood of the inward crossing 
where the classical geometry is a good approximation can be large. It seems that 
sufficiently large time delays would allow for arbitrarily slow change of colour. We 
cannot exclude, therefore, that the quantum spacetime contains an extended region 
with the geometry resembling its classical counterpart near a black hole horizon, 
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at least locally. This can be true even if the quantum spacetime as a whole differs 
strongly from any typical classical collapse solution. 

One can even imagine the following scenario (which needs a more realistic model 
than a single thin shell). A quantum system with a large energy collapses and 
re-expands with huge time delay. The black hole horizon phase is so long, that 
Hawking evaporation becomes significant and must be taken into account in the 
calculation. It does then influence the time delay and the period of validity of the 
black hole approximation. The black hole becomes very small and only then the 
change of horizon colour becomes significant. The white hole stage is quite short 
and it is only the small remnant of the system that, finally, re-expands. The whole 
process can still preserve unitarity. In fact, this is a scenario for the issue of Hawking 
evaporation process. At least, it is not excluded by the results of the present paper. 

The calculations of this paper are valid only for null shells. Similar calculations 
have been performed in f2(| . There has been re-expansion and unitarity for massive 
shells if the rest mass has been smaller than the Planck mass (10~5 g). It is very 
plausible that the interpretation of these results is similar to that given in the present 
paper. Thus, we can expect the results valid at least for all "light" shells. There is, 
in any case, a long way to any astrophysically significant system and a lot of work 
is to be done before we can claim some understanding of the collapse problem. 

Our method of dealing with the problem employs simplified models and a kind 
of effective theory of gravity; it does not worry about the final form of a full-fledged 
theory of quantum gravity. This need not be completely unreasonable approach. 
Even if the ultimate quantum gravity theory were known, most calculations would 
still be performed within the approximation of some effective theory and for simpli- 
fied models (compare the situation in the QCD). The method can give useful hints 
also because of the fact that the black hole geometry is "made up" from purely de- 
pendent degrees of freedom of the gravitational field, and these degrees of freedom 
have no proper quantum character of their own. 

To summarize: We have demonstrated, at least for light shells, that quantum 
theory can smoothly unify two states of motion, one being the time reversal of the 
other, into one history. In this way, geometry containing a piece of a black hole 
horizon can be followed by geometry containing a piece of a white hole horizon — 
just the opposite to the situation we know from the Kruskal diagram of the classical 
general relativity. In this way, the quantum evolution can stay unitary and the 
question posed at the beginning of the paper can be answered as follows: A quantum 
system is not always lost if it falls under its Schwarzschild radius. 
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